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Anisotropic lattices are an efficient means of studying the glueballs of QCD, however problems arise with 
simulations of the lightest, scalar state. The mass is strongly dependent on the lattice spacing, even when a 
mean-held improved gluon action is used. The nature and cause of these errors are discussed and the scaling 
properties of the scalar from different lattice actions are presented. 


1. THE SCALAR DIP 

Efficient studies of the QCD glueball spectrum 
are made possible on coarse spatial lattices, pro¬ 
vided the temporal discretisation is finer. This 
permits sufficient measurements of the glueball 
correlator before statistical fluctuations dominate 
(^. These studies found good scaling for the 
2++ and I"' glueballs. However for the lightest, 
scalar state, the mass (in units of ro) was seen to 
fall rapidly to a minimum when the spatial lattice 
spacing is 0.25 fm, where the scaling violations 
are « 25% and then rise as the lattice spacing 
is increased further; the “scalar dip”. This effect 
was also seen in SU{2) simulations [H, although 
the dip is less pronounced. 

The first question is whether Symanzik im¬ 
provement of the action has worked. Fig. |^shows 
the scaling of the scalar from both Wilson and 
improved simulations. In this figure, the mass 
of the U-irrep tensor glueball is used to set the 
scale as these mass ratios are independent of any 
anisotropy renormalisations. The two curves are 
results of fitting the datasets to a parabola. Fig. 
demonstrates the improvement reduces the depth 
of the dip, but only by about a third. 

We are examining a number of other explana¬ 
tions for this behaviour; the first is that using 
the plaquette definition for the mean-link param¬ 
eters leads to an under-estimate of the effects of 
the tadpole graphs. Secondly, the influence of a 
critical end-point in a higher-dimensional space 
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of couplings may be at fault. 

In our preliminary study of the mean-link def¬ 
inition no significant improvement is observed. 
This is in contrast to findings for SU{2) where 
the weaker scalar dip was reduced further by em¬ 
ploying the Landau-link definition Q. 

Simulations of QCD in the fundamental-adjoint 
space of plaquette couplings demonstrate the ex¬ 
istence of a critical point close to the fundamental 
axis. As this point is approached, the mass of the 
scalar (in units of the string tension) has been 
shown to fall rapidly j|]. At the critical point, at 
least one lattice correlation length must diverge 
and this is seen to be the scalar channel. If this ef¬ 
fect is responsible, addition of negatively coupled 
adjoint terms to the action might define a scal¬ 
ing trajectory where the influence of this critical 
point is diminished Q. 

This solution ought to be independent of 
any Symanzik-improvement programme. To test 
this, we have added adjoint-like terms to the 
anisotropic Wilson action. This report dis¬ 
cusses progress with these studies. Simulations of 
Symanzik-improved actions with these new terms 
are in progress. 

2. DESIGNING NEW ACTIONS 

Based on the observation of the nearby critical 
point, we consider one member of a class of new 
actions which include “adjoint-like” terms. The 
term we add is designed to be easy to simulate 
with existing update methods. 




Figure 1. The effect of Symanzik improvement 
on the scalar glueball. 


Beginning with the plaquette operator, 

P^Li^ix) = ^R-eTr Uij^{x)U^{x+ij,)Ul{x+i>)Ul{x) 

the Wilson discretisations of the magnetic and 
electric field strengths are 

ns=J2{l- Pij{x)} = I /a Tr + 0{a^) 
and 

n* = V {1 - P*o(x)} = ^ /a Tr £2 + 0{a^) 

respectively, where i,j are spatial indices and ^ is 
the tree-level anisotropy, ag/at- We introduce a 
term which correlates pairs of spatial plaquettes 
separated by one site temporally 

= ^ E {i-a,(x)p.,(x + t)}. 

x4>j 

A similar term with spatially correlated plaquette 
pairs can be constructed but is not considered 
here. 

The separation of the two plaquettes allows the 
standard Cabibbo-Marinari and over-relaxation 
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Figure 2. Glueball mass ratio dependence on 
the two-plaquette mixing parameter u for the 
and irreps at fixed p. The 

horizontal bands indicate the continuum limit. 
Dashed lines guide the eye. 


gauge-field update methods to be applied. In¬ 
cluding two-plaquette terms adds a computa¬ 
tional overhead of just 10% to our anisotropic 
Wilson action workstation codes. 

It can be shown that the operator combination 

fig = (1 -I- w) fig — w 

has an identical expansion in powers of Og^t to fig 
up to 0(ag j). Thus a signal for being closer to 
the QCD fixed point would be physical ratios be¬ 
coming weakly dependent on the free parameter, 
to. 

A candidate action (without Symanzik im¬ 
provement) is then 

i ■ 

Mean-link improvement is applied in the standard 
manner. 

2.1. Simulation results 

We have performed a range of simulations at a 
number of values of w and lattice spacings. 







Figure 3. Scaling properties of the scalar of Sui- 





Figure 4. Scaling properties of the and ^ 
irreps. Filled symbols indicate data. 


Fig. H shows the dependence of ratios of glue- 
ball masses on oj at fixed /3 (= 2.7) for Si^. Close 
to tu = 0 (anisotropic Wilson), these scaling ra¬ 
tios show a strong dependence on w, both in the 
scalar and pseudoscalar channels. In the range 
a; = 1 to 3, the dependence is much weaker for all 
three channels, and the ratios are closer to their 
continuum values. The lattice spacing does not 
change significantly over this range of w; the lat¬ 
tice mass of the E irrep varies by just 7%. 

Fig. ^compares the scaling properties of S^j for 
two values of uj to the anisotropic Wilson action. 
The cut-off dependence in the scalar channel is 
seen to be reduced significantly. More data are 
required to resolve the continuum limit fully. 

Fig. ^ shows the tensor irrep rotational symme¬ 
try violations and the scaling of the pseudoscalar 
glueball. Remarkably, the large scaling violations 
of the Wilson action in these channels are seen al¬ 
most to vanish for with w = 1 or 2. This was 
not anticipated from critical point arguments. 

3. CONCLUSIONS 

By comparing data from anisotropic Wilson 
simulations to their Symanzik-improved counter¬ 


parts, we have shown that improvement reduces 
the scaling violations of the scalar glueball, how¬ 
ever large effects remain. 

We found that adding to the action a new 
term which includes the product of two plaquette 
traces leads to significant reduction in the curi¬ 
ous finite-lattice-spacing effects seen for the scalar 
glueball mass. Surprisingly, the pseudoscalar and 
T 2 tensor irrep masses (in units of the E-irrep 
mass) also show far less cut-off dependence. We 
emphasise that the action, S^j used here is not 
Symanzik improved, and has leading discretisa¬ 
tion errors at 0{al^). 

We are performing simulations where two- 
plaquette terms are added to Symanzik-improved 
actions of the type considered in Ref. . 
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